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Abstract 

We consider dynamics of Josephson vortex lattice in layered superconductors with magnetic, 
charge (electrostatic) and charge-imbalance (quasiparticle) interactions between interlayer Joseph- 
son junctions taken into account. The macroscopic dynamical equations for interlayer Josephson 
phase differences, intralayer charge and electron-hole imbalance are obtained and used for numeri- 
cal simulations. Different transformations of the vortex lattice structure are observed. It is shown 
that the additional dissipation due to the charge imbalance relaxation leads to the stability of 
triangular lattice. 
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The investigation of the resistive state of layered high-T c superconductors in the external 
magnetic field parallel to the layers is the important task of modern experiment and theory. 
The c-axis resistivity in this case is determined by Josephson flux flow. Thus, the nonsta- 
tionary properties of HTSC can be investigated in wide parameter range of temperatures, 
fields and currents. In the recent experiments |IJ, [3], |J flux flow states in low magnetic 
fields are observed and in Ref. []5], - in high fields up to 3.5 T. The number of peculiarities is 
discovered. In particular, in Ref. || [| multiple flux flow branches, associated with different 
vortex modes, are observed. In Ref. @, || it is shown that two universal flux flow regimes 
with different V/ H take place at different magnetic fields and transformation from the first 
state to the second is possible with current rise at intermediate fields. Finally, in Ref. 
broad-band microwave emission is measured with frequencies much less than Josephson fre- 
quency ujj = 2eV/Nh. All this phenomena are related to the dynamical transformations of 
Josephson vortex lattice (JVL) and vortex interaction with linear plasma modes (see also 

Ref. 0, 0, H, 01 0)- 

On the other hand, there is the long-standing problem to obtain the coherent electromag- 
netic emission from stacked Josephson junction. The most simple way to do this is to achieve 
"in-phase" (square) arrangement of vortices. Although at zero external current the triangu- 
lar JVL is favorable, it was shown for two-junction stack theoretically and experimentally 
n| , |TJ|] and for multi-junction stack theoretically || [II], O, [T5|, |TB| that moving square JVL 



can be stable. However no indications of this in-phase regime were found in experiments 
with HTSC. The reason for this discrepancy may be that only magnetic coupling between 
layers have been taken into account in calculations. 

In [0, 0, |(], |2l], H|, [23j it was also suggested that in the case of thin superconduct- 
ing layers some nonequilibrium mechanisms are to be taken into account, namely electrical 
charging of layers (charge effect) and nonequilibrium quasiparticles (charge-imbalance ef- 
fect). To investigate this effects we obtain macroscopic dynamical equations for interlayer 
Josephson phase differences with magnetic, charge and charge-imbalance interactions taken 
into account. 

Following |7], [8| we use London expression for the longitudinal in-layer supercur- 
rent neglecting longitudinal in-layer electric field and normal current (with accuracy 
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together with the equation B = rotA and the definition of invariant Josephson phase 
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here B is the averaged interlayer magnetic field, z is the unity vector perpendicular to the 
layer and d is the interlayer distance, layers assumed to be thin (d <C 



Then, using Maxwell equation rotB = + ^fj 3 projected on to z-axes, we obtain 
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Taking use of the continuity equation 2fia + Viijj[ + in '" +1 in ' 1 



edivE = 47rp in integral form p 



()/ v ||,,„ ( , — and Maxwell equation 
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where j ex t is the external current in overlapping geometry (for the details of introducing the 
external current in different geometries see Ref. |7|]). 

To find the interlayer Josephson current j n n +i we use the theory of nonequilibrium 
Josephson effect developed recently (|Tj| [L^, |2(| [21], and references therein). We 

take into account that in nonequilibrium state there is nonzero invariant potential 



$ n {t) = (f> n + {h/2e){d6 n /dt), 
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where n is the electrostatic potential and 8 n is the phase of superconducting condensate, 
$ = in the equilibrium state. 

An ordinary Josephson relation (dip/dt) = (2e/H)V between the Josephson phase differ- 
ence t£ n ,n+i and voltage V n , n+ i = E n>n+ xd is violated. Instead, we have (from (|5])) 
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The charge density inside a superconducting layer is 

p n = -2e 2 N (<5> n - tf n ), 



(8) 



where the first term is the charge of superconducting electrons determined by the shift of the 
condensate chemical potential 5p s = — e$ and the second term is the quasiparticle charge, 
described by potential (details can be found in Ref. J23J). Finally, the equation for the 
charge imbalance can be written in the form corresponding to the generalized nonstationary 
Ginzburg-Landau theory ]23f 



^ = (1 - T)— + 2v t - — j + 2i*(* 4 _ 1 + % +1 - 2%) - r q (9) 

where rj = 2v t r q is the parameter of disequilibrium, r q is the well-known charge-imbalance 
relaxation time, v t = (Ae^NoRSdo) -1 is the "tunnel frequency" , R is the normal resistivity of 
the tunnel junction, V = Sdo is the volume of the superconducting layer, N = mp^/27r 2 /l 3 
is the density of state (or, instead, Nf = N do can be considered as 2D density of states), 
T <C 1 is the parameter dependent on the details of microscopic theory. 

In the same approximation we derive the expression for the nonequilibrium interlayer 
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JM+i = Jc sm <p iji+1 + + — r— (10) 

Finally, in dimensionless form we obtain 

Vf = j* >i+1 ~ S(j*_ hi + j* +14+2 ), (11) 
(9Vi,i+l . dlf^i+x ,a( d Vi 9Pi+l\ ■ / 10 x 

= ^ g r 2 + -^T" + sm Wi+i + ^ - ^+1 + Z 3 I 7^7 - -^T" J - (12) 
a ^ + ^i + ViWi ~ ^i-i - ^»+i) = ^ I — ^1 I + a (1 - r )^T' ( 13 ) 

A*i + C(2/ii - - = % + C I — ^- I • (I 4 ) 

where //(t) = $(t)/V c , ^(t) = #(*)/K, K = fc<j c /2e, a' = r^ c , /5 = (fteu; 2 S)/(87re$ c ), 
Aj = (/ic 2 do)/(167rej c (A 2 + dd )), ( = (erj)/(d d), uj c = 2eRj c /h, r = u c t, x = x/Xj, 
s = 0.5/(1 + (dd )/(2X 2 ))). 

This system of equations is a generalization of well-known magnetic coupling model of 
layered superconductors [|7|, ||. Neglecting electron-hole imbalance (77 = 0, \l/ n = 0), we 
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derive the system of equations obtained in Ref. fL5f , and neglecting the charge effect (( = 0, 
\l/ n = <£> n ) we derive the system of equations similar to that of |18], [19], plf . 

Before further study, the Eqs. (pTT])- (|14[) must be added with boundary conditions, which 
we write as 

J -1,0 = Jn,n+i = 0) 
= B, 



dip 



3X 0,L 
IpO = = 



(15) 



dx 



0. 



0,L 

Ak> = = 0. 

For numerical solution of the equations (1 1 1)-(14) with boundary conditions ([TB]) we in- 
troduce the new functions <p, v, z, defined by the relations 



0i{x,t) = ip i;i+1 (x,t) - Bx, 

Vi (X, t) = -^Tj^ - + & fa, t) ~ (X, t) , 

Zi(x,t) = £ ^fl^L _|_ a '(i)i(x,t) - ijj i+l (x,t)), 
i = 0,...N-1, e = a'(l-r). 
Let us denote now by U the column vector-function of kind 
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(here is the difference Laplasian, (Adf)i = fi-i — 2 fa + fi+i) and by F the column vector 
of kind 
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Then we can rewrite the above system of equations in the following "operator" form: 

dt (20) 
(sA d + 2s- l)f = -Vjj& 

with appropriate homogeneous boundary conditions. For numerical solution of the equation 
( p0|) we use, after sampling in x-direction, the semi-implicit scheme of form 

U(t + At) = 2(E - ^A)-\U(t) + jF(t)) - U(t) - jF(t), (21) 

in which all needed inverse operators can be easily evaluated by standard sweep method. 
The transformation from vector U to initial variables is evident. 

The developed program gives us the possibility to observe the state of system during 
calculation process in convenient graphical form. We can therefore find many different 
variants of system dynamics and vortex lattice structures. 

The most part of simulations has been carried out for the case of high magnetic fields 
H > H* = $o/7^ 2 , where $o = irfrc/e is the flux quantum, 7 = \ c /\ a b is the anisotropy 
parameter, and t = d + d is the period of the structure. At this fields triangular vortex 
lattice is formed in the static case. 

The results of numerical simulations can be shortly described in the following way. At 
small currents vortex lattice remains triangular. At high currents the situation depends on 
dissipation. If interlayer dissipation is strong enough (/? < 1) then triangular lattice is stable 
at all currents. Otherwise, we observed various transformations of lattice structure (further 
(3 = 100 and s = 0.48 are considered). 

(i) At pure magnetic coupling (£ = 0, r\ = 0) we obtain the typical picture of JVL 
transformations. The example is shown on Fig. 1. We consecutively change current from 
zero to some high value {a-f) and then decrease it up to zero again (g-l). The square lattice 
is clear seen on Fig. d,e,i, and j. Low current regimes a,l, as well as high current one f,g are 
triangular (or close to it). Regimes b,c and k are " inhomogeneous" , - velocities of vortex 
chains are different in different layers. These regimes are also "breathing" due to strongly 
excited nonlinear plasma waves. We find that such breathing modes can be triangular or 
square in average. Note that different regimes exist at the same current that can lead to 
hysteresis on VCC. 

(ii) If we take charge coupling into account (Fig. 2, ( = 1, r] = 0), the picture of JVL 
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transformations is qualitatively the same but the regimes with square lattice are shifted to 
higher currents. 

(iii) For charge-imbalance (quasiparticle) coupling the result is qualitatively different. We 
have found that if the time of charge-imbalance relaxation is large enough (approximately 
a' = T q uj c > 1) than the triangular lattice become stable even if the parameter of disequi- 
librium rj is small. The example is shown on Fig. 3 (( = 0.1, rj = 0.1, a' = 100). There 
is no JVL transformations at this parameters although both charge and charge-imbalance 
couplings are weak. The origin of this effect is the additional dissipation due to the charge- 
imbalance relaxation. The detailed theory will be considered in other publication. 

In conclusion, we derive macroscopic dynamical equations for layered superconductors 
with magnetic, charge and nonequilibrium quasiparticle interactions taken into account. 
Many dynamical regimes are observed, in particular triangular, square, inhomogeneous and 
breathing modes. We established that additional dissipation due to the charge-imbalance 
relaxation can prevent JVL transformations and make triangular lattice stable at all currents. 
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Figure captions 



Fig.l. The structure of Josephson vortex lattice (sin</? n are shown) with (3 — 100, s — 
0.48, rj = 0, C = 0, j = 0.2, 0.4, 0.6, 0.8, 1, 1.2. Current is increased (a-f) and then decreased 
(9-1)- 

Fig.2. The structure of Josephson vortex lattice with j3 = 100, s = 0.48, rj — 0, ( — 1. 
Fig.3. The structure of Josephson vortex lattice with /3 = 100, s = 0.48, rj = 0.1, C = 0.1, 
a' = 100, T = 0.01. 
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